Determination of the coherence length and the Cooper-pair size in unconventional 

superconductors by tunnelling spectroscopy 
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The main purpose of the paper is to discuss a possibility of the determination of the values of the 
coherence length and the Cooper-pair size in unconventional superconductors by using tunnelling 
spectroscopy. In the mixed state of type-II superconductors, an applied magnetic field penetrates 
the superconductor in the form of vortices which form a regular lattice. In unconventional super- 
conductors, the inner structure of a vortex core has a complex structure which is determined by the 
order parameter of the superconducting state and by the pairing wavefunction of the Cooper pairs. 
In clean superconductors, the spatial variations of the order parameter and the pairing wavefunction 
occur over the distances of the order of the coherence length and the Cooper-pair size, respectively. 
Therefore, by performing tunnelling spectroscopy along a line passing through a vortex core, one is 
able, in principle, to estimate the values of the coherent length and the Cooper-pair size. 

PACS numbers: 74.25.Op; 74.25. Qt; 74.50. +r; 74.25.Jb 



I. INTRODUCTION 

Superconductivity requires the electron pairing and 
the onset of long-range phase coherence. These two phys- 
ical phenomena are independent of one another. In the 
framework of the Ginzburg-Landau theory [1], the su- 
perconducting state is characterized by the order param- 
eter 'J'. The coherence length £cl is the characteris- 
tic scale over which variations of '5 occur, for example, 
near a superconductor-normal metal boundary. Gener- 
ally speaking, the coherence length is different from the 
size of the Cooper pairs, £, which is related to the wave- 
function of a Cooper pair ip. Furthermore, the coher- 
ence length depends on temperature, £gl(T), while the 
Cooper-pair size is temperature- independent, at least, in 
conventional superconductors. The coherence length di- 
verges at T — > T c , where T c is the critical temperature. 
For every superconductor, the knowledge of the values 
of the coherence length and the size of Copper pairs is 
important for the understanding of the underlying mech- 
anism of superconductivity. 

The superconducting state can be destroyed by a suf- 
ficiently strong magnetic field. The variation of the ther- 
modynamic critical field H c with temperature for a type- 
I superconductor is approximately parabolic: H C (T) ~ 
JT C (0)[1 - (T/T c ) 2 ], where H c {0) is the value of the crit- 
ical field at absolute zero. For a type-II superconductor, 
there are two critical fields, the lower critical field H c \ 
and the upper critical field H C 2 . In applied fields less than 
H c \, the superconductor completely expels the field, just 
as a type-I superconductor does below H c . At fields just 
above H c i, flux, however, begins to penetrate the super- 
conductor in microscopic filaments called vortices which 
form a regular lattice, as shown in Fig. 1(a). Each vortex 
consists of a normal core in which the magnetic field is 
large, surrounded by a superconducting region, and can 



be approximated by a long cylinder with its axis parallel 
to the external magnetic field. Inside the cylinder, the 
superconducting order parameter is zero, as illustrated 
in Fig. 1(b). The radius of the cylinder is of the order 
of the coherence length £gl- The supercurrent circulates 
around the vortex within an area of radius ~ A, the pen- 
etration depth. The vortex state of a superconductor, 
predicted theoretically by Abrikosov [2], is also known 
as the mixed state. By increasing the magnitude of the 
applied magnetic field from H c \ to H C 2, the distance be- 
tween vortices decreases, becoming zero at H c i- At H C 2, 
the field penetrates completely the superconductor, mak- 
ing it normal. 

Probably, the best technique for the determination of 
the behaviour of the vortex cores in the bulk of type-II 
superconductors is muon spin rotation (/xSR) measure- 
ments [3]. However, /xSR are not able to determine (i) 
the inner structure of the vortex cores and (ii) the local 
density of states of quasiparticle excitations. Generally 
speaking, there is no perfect technique: every technique 
has some disadvantages. 

Tunnelling spectroscopy is an unique probe of the su- 
perconducting state in that it can, in principle, reveal 
the quasiparticle excitation density of states directly with 
high energy resolution [4] . By using a scanning tunnelling 
microscope (STM) , one can obtain images of the vortex 
lattice in the mixed state of superconductors, and per- 
form local spectroscopy inside and outside vortex cores 
[5-14]. At low magnetic field H ~ H c \, when the dis- 
tance between vortices is large (3> £gl), the variation of 
the local density of states inside and outside vortex cores 
is determined by the order parameter '5 and the pairing 
wavefunction of the Cooper pairs ip. The spatial varia- 
tions of \& and ip occur over the distances of the order 
of £gz, an d £, respectively (in the clean limit). There- 
fore, by performing tunnelling spectroscopy along a line 
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FIG. 1. (a) In the mixed state of type-II superconductors, 
an applied magnetic field penetrate the superconductor in the 
form of vortices which form a regular lattice. Each vortex 
consists of a normal core, and can be approximated by a long 
cylinder with its axis parallel to the external magnetic field. 
The radius of the cylinder is of the order of the coherence 
length £g£- (b) Spatial variations of the magnetic field H and 
the order parameter ^ inside and outside an isolated vortex in 
an infinite superconductor. R is the distance from the center 
of the vortex, and £gl and A are the coherence length and 
the penetration depth of the superconductor, respectively (in 
type-II superconductors, £gl < A). 



passing through the center of a vortex core, one can in 
principle estimate the values of the coherent length and 
the Cooper-pair size. The main purpose of the paper is 
to investigate a possibility of the determination of the 
values of the coherence length and the Cooper-pair size 
in unconventional superconductors by using tunnelling 
spectroscopy in the mixed state. 

The paper is organized as follows. We first discuss 
the spatial variation of the local density of states in the 
vortex cores of conventional superconductors. Then, we 
shall turn our attention to half-conventional and uncon- 
ventional superconductors. In the following section, we 
shall compare our results with real data obtained in half- 
conventional and unconventional superconductors. The 
paper ends with a discussion and conclusions. 



II. CONVENTIONAL SUPERCONDUCTORS 

Superconductivity requires the electron pairing and 
the onset of long-range phase coherence. In the frame- 
work of the BCS (Bardeen-Cooper-Schrieffer) theory for 
conventional superconductors [15], the electrons form 
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FIG. 2. Temperature dependences of (a) the coherence 
length £gl and the Cooper-pair size £, and (b) 1/£gl and l/£ 
in "clean" conventional superconductors. In a first approxi- 
mation, the Cooper-pair size is independent of temperature. 
The dependences £gl(T) and 1/£gl(T) are shown schemati- 
cally. 

pairs due to phonons, while the phase coherence is es- 
tablished by the overlap of the Cooper-pair wavefunc- 
tions. The latter process is also called the Josephson 
coupling. Such a way of the establishment of the long- 
range phase coherence gives rise to an order parameter 
which is a "magnified" version of the Cooper-pair wave- 
functions. Therefore, the values of the coherence length 
and the Cooper-pair size in conventional superconduc- 
tors coincide at T — 0: £gl(0) = £(0) = as shown 
in Fig. 2(a). The coherence length (or the Cooper-pair 
size) £o determined by the energy gap at zero tempera- 
ture, A(T = 0), is called intrinsic: 
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(1) 



where vf is the Fermi velocity (on the Fermi surface), and 
h = h/2ir is the Planck constant. Let us estimate £o- I n a 
metal superconductor, A(0) ~ 1 meV. Substituting this 
value into Eq. (1), together with vp ~ 1.5 x 10 8 cm/s and 
h = h/2n ~ 6.5 x 1(T 13 meVs, we obtain ( ~3x 1CT 5 
cm = 3xl0 3 A. Thus, the values of the intrinsic coherence 
length in conventional superconductors is very large in 
comparison with the interatomic distances (~ 1-2 A). In 
conventional superconductors, the electron pairing and 
the onset of phase coherence occur simultaneously at T c . 
It is worth mention that Eq. (1) can be obtained in a 
good approximation by using the uncertainty principle 
Sp^o ~ h, where Sp is the momentum spread [16]. In 
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conventional superconductors, the assumption that £(T) 
is constant with temperature is sufficiently good, because 
the T c value in conventional superconductors is very low 
and there are no structural transitions below T c . 

Consider the so-called clean and dirty limits for super- 
conductors. In the clean limit, I ^ £o> where I is the 
electron mean free path; in the dirty limit, I <C £o- At 
< T < T c , the coherence length in "clean" conven- 
tional superconductors is always larger than the size of 
Cooper pairs, £ < £gl(T), as depicted in Fig. 2(a). In 
the framework of the Ginzburg-Landau theory, the tem- 
perature dependence of the coherence length in "clean" 
superconductors at temperatures close to T c is given by 



& L {T)=0.74£ 
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T 
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(2) 



For "dirty" superconductors, the Ginzburg-Landau tem- 
perature dependence of the coherence length at T ~ T c 



^l(T) = 0.85(^) 1/2 



T 



-1/2 



(3) 



From Eqs. (2) and (3), one can see that £g L — > oo as 
T — > T c . It is worth to mention that, in very "dirty" 
metals, the mean electron free path £ plays the role of 
the coherence length [17]. 

The condensation of the Cooper pairs in conventional 
superconductors, and therefore, of the condensate, oc- 
curs in momentum space. In this case, one can argue that 
the "real" characteristics of the superconducting conden- 
sate in conventional superconductors are not £gl and £ 
but and £ _1 shown in Fig. 2(b). The temperature 
dependence of ^q 1 l in Fig. 2(b) is similar to the BCS 
temperature dependence of the energy gap, A(T), illus- 
trating the increase of the phase stiffness with decreasing 
temperature. 

Let us now discuss the variations of the characteris- 
tics of the superconducting state into the vortex cores 
appearing in the mixed phase of conventional supercon- 
ductors. The spatial variations of the magnetic field and 
the order parameter inside and outside an isolated vortex 
are illustrated in Fig. 1(b). As shown by Gor'kov, the 
spatial variations of the energy gap in conventional super- 
conductors are proportional to the variations of the order 
parameter, A cx |^| [18]. Therefore, in conventional su- 
perconductors, the energy gap also goes to zero inside the 
vortex core, as shown in Fig. 3. In conventional super- 
conductors, by performing tunnelling spectroscopy along 
a line passing through a vortex core, it is possible to esti- 
mate only the value of the coherence length £gl but not 
£. On the other hand, in conventional superconductors, 
£ ~ £ Gi at T < T c /2. As a consequence, the knowledge 
of ^gl at low temperatures is sufficient. 

Alternatively, the value of the coherence length can be 
estimated in momentum space. By performing topogra- 
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FIG. 3. Spatial variations of the energy gap A inside and 
outside an isolated vortex in conventional superconductors. 
R is the distance from the center of a vortex core. 

phy image above vortex cores at bias |V| < A/e, where 
e is the electron charge, and then, making the Fourier 
transform of the image, one can in principle estimate the 
value of the coherence length in accordance with the plot 
in Fig. 2(b). 

In a "clean" conventional, s-wave superconductor, 
quasiparticles whose energy E is less than the bulk en- 
ergy gap, E < A, may theoretically form discrete local- 
ized states in the vortex core [19,20]. Andreev reflections 
of quasiparticles from the normal core-superconductor in- 
terface give rise to these bound states. The lowest energy 
of the bound states is approximately ~ A 2 /2Ep, where 
Ep is the Fermi energy [19]. In a superconductor with 
an order parameter having nodes, the bound state will 
"leak out" through the nodes, giving rise a broad peak at 
the Fermi level in the quasiparticle local density of states 
of vortex cores [21]. 



III. HALF-CONVENTIONAL 
SUPERCONDUCTORS 

In this section we discuss the mixed phase in two- 
band superconductors. The Fermi surface of such su- 
perconductors consists of, at least, two disconnected sec- 
tions. Below a certain temperature, the Cooper pairs 
are formed in one section of the Fermi surface due to 
the electron-phonon interaction. Then, due to either 
the interband scattering or Cooper-pair tunnelling, or 
both, the electron pairing is induced into the second 
section (band) of the Fermi surface. So, there are two 
types of the Cooper pairs in the system, and all of 
them are formed due to the electron-phonon interaction. 
Therefore, we shall call such superconductors as half- 
conventional. Consider a system in which the genuine 
Cooper pairs formed in the first band cannot establish 
alone the long-phase coherence. This can happen in a 
system in which the size of the Cooper pairs is small (in 
comparison with that of the Cooper pairs in metals) and 
their concentration is low. As a consequence, their wave- 
functions do not overlap steadily. The long-range phase 
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FIG. 4. Temperature dependences of the coherence length 
£gl and the Cooper-pair size £l and £s in "clean" two-band 
superconductors. In two-band superconductors, there two 
types of the Cooper pairs having the different size. In a first 
approximation, £l and £s are independent of temperature. 
The dependence £gl(T) is shown schematically. 

coherence of the superconducting state appears owing to 
the Josephson coupling of the wavefunctions of the in- 
duced Cooper pairs. In reality, in two-band superconduc- 
tors the genuine Cooper pairs are low-dimensional, i.e. 
one- or two-dimensional. At the same time, the induced 
Cooper pairs in the second band are three-dimensional. 

Let us label the energy gap of the low-dimensional 
Cooper pairs by A^, and their size by £l. Analogously, 
for the three-dimensional Cooper pairs, Ag and £s- The 
labels "L" and "S" stand for "Large" and "Small", re- 
spectively, because always A^ > As- This is typical for 
induced superconductivity [17]. In contrast, £l < £g, 
i.e. the size of the three-dimensional Cooper pairs is al- 
ways larger than that of the low-dimensional ones. This 
follows from the fact that £ oc 1/ A for any type of super- 
conductivity. Denote the wavefunctions of the two types 
of the pairs by ipi, and tps- The order parameter of the 
superconducting state \& in the system is a "magnified" 
version of ipg, and the coherence length ^ql of the order 
of £5, as schematically shown in Fig. 4. 

We are going now to consider the variations of the 
characteristics of the superconducting state into the vor- 
tex cores appearing in the mixed state of "clean" half- 
conventional superconductors. Figure 5(a) shows the 
spatial variations of the order parameter and the aver- 
aged wavefunction into a vortex core. Inside a vor- 
tex core, the order parameter always goes to zero, while 
\ip L \ should not. Depending on the coupling strength of 
the low-dimensional Cooper pairs, there are three possi- 
ble outcomes for the behaviour of \ into a vortex core, 
shown by the letters A, B and C in Fig. 5(a). The spa- 
tial variation of the averaged wavefunction, if such 
occurs, is of the order of £l- From the discussion in the 
previous section, the smaller energy gap Ag will be al- 
ways zero inside a vortex core, as depicted in Fig. 5(b). 
In contrast, the larger energy gap A^ should not be zero 
into a vortex core, and even, may remain unchanged [the 
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FIG. 5. Spatial variations of (a) the order parameter |^| 2 
and the averaged wavefunction IV'il 2 , and (b) the energy gaps 
Al and As inside and outside an isolated vortex in an infinite 
"clean" two-band superconductor. R is the distance from the 
center of the vortex. £gl is the coherence length, and £l is 
the size of the smaller (genuine) Cooper pairs (£l < £gl). 
For simplicity, all the functions are normalized at R = 00 to 
1. Since the pairing wavefunction ipL is not global, ipL in plot 
(a) is average. The letters A, B and C indicate a different 
behaviour of IV'lI 2 and Al. 

case A in Fig. 5(b)]. In the latter case, the size of the 
Cooper pairs cannot be determined from the tunnelling 
data. Therefore, we discuss further only the cases B and 
C sketched in Fig. 5(b). The case C in Fig. 5 can theoret- 
ically occur; however, it is less likely that such a situation 
can be realized in practice because the existence of the 
induced Cooper pairs depends on the presence of the gen- 
uine Cooper pairs. In vortex cores of half-conventional 
superconductors, tunnelling measurements tuned to As 
should show the disappearance of the energy gap in the 
center of a vortex core, as illustrated in Fig. 6(a). Prob- 
ing Ai, one should observe its variation inside a vortex 
core sketched in Fig. 6(b). Since this variation of Al oc- 
curs over the distance of the order of £l , one can estimate 
the value of £l from the tunnelling data. 

It is important to note that, in two-band superconduc- 
tors in which the genuine Cooper pairs can along estab- 
lish the long-range phase coherence, the magnitudes of 
both Ag and A^ will be zero in the center of a vortex 
core. Such a situation was discussed elsewhere [22]. To 



4 



--vortex core-- 









(a) 










A L 










(b) 






V c / 








R 



FIG. 6. Spatial variations of the energy gaps (a) As and 
(b) Ai inside and outside an isolated vortex in an infinite 
"clean" two-band superconductor. The line R passes through 
the center of the vortex core. £gl is the coherence length, 
and £t is the size of the smaller (genuine) Cooper pairs 
(£l < £gl). For Al, the letters B and C correspond to the 
cases shown in Fig. 5(b). 

remind, we consider here two-band superconductors in 
which the genuine Cooper pairs formed in the first band 
cannot establish alone the long-phase coherence. The 
long-range phase coherence in this type of superconduc- 
tors appears due to the overlap of the wavefunctions of 
the induced Cooper pairs in the second band. 

As in the case of conventional superconductors, in half- 
conventional superconductors it is also not possible to 
estimate the value of £s from the tunnelling spectroscopy 
data. On the other hand, £s ~ £gl at T < T c /2. For the 
estimation of the value of £l, a hint can be suggested. It 
can be estimated from the following equation 

&~£gl^. (4) 

In half-conventional superconductors, the two energy 
gaps have both an s-wave symmetry typical for the 
electron pairing due to phonons. They are both 
anisotropic: if A5 is only slightly anisotropic, Al is 
highly anisotropic, and can even have nodes. In half- 
conventional superconductors, the maximum magnitude 
of A^ is a few times (~ 3) larger than that of Ag. This 
is typical for induced superconductivity [17]. 

IV. UNCONVENTIONAL SUPERCONDUCTORS 

Before we discuss the mixed phase of unconventional 
superconductors, let us first give the definition of an 
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FIG. 7. Temperature dependences of the coherence length 
£gl and the Cooper-pair size £ in unconventional supercon- 
ductors in which the superconducting condensate has a struc- 
ture similar to that of a Russian doll (see text for more de- 
tails). T p is the pairing temperature. In plot, £ is independent 
of temperature; however, in some cases, this assumption can 
be incorrect. The dependence £gl{T) is shown schematically. 

unconventional superconductor. In a wide sense, a su- 
perconductor is unconventional if the mechanism of su- 
perconductivity in this superconductor is different from 
the BCS mechanism and from that of half-conventional 
superconductors. In this section, we shall consider two 
types of unconventional superconductors. We start with 
superconductors in which the superconducting conden- 
sate has a structure similar to that of a Russian doll. 

A. Russian-doll-like condensate 

In these superconductors, both the quasiparticle pair- 
ing and the long-range phase coherence are mediated by 
bosonic excitations. If the electron pairing is always me- 
diated by bosonic excitations, this is not typical for the 
establishment of the phase coherence. In conventional 
and half-conventional superconductors, the long-range 
phase coherence occurs due to the overlap of the Cooper- 
pair wavefunctions, thus, without participation of any 
bosonic excitations. In superconductors with a Russian- 
doll-like condensate, the electron pairing is mediated, for 
example, by phonons, and the phase coherence, for in- 
stance, occurs due to spin fluctuations. Assume that the 
Cooper pairs are formed above the critical temperature, 
T p > T c . In this type of unconventional superconduc- 
tors, the symmetry of the Cooper-pair wavefunctions can 
differ from the symmetry of the order parameter of the 
superconducting condensate. Thus, the superconducting 
condensate in this type of superconductors has a struc- 
ture similar to that of a Russian doll. 

Figure 7 shows the temperature dependences of the co- 
herence length £gl and the Cooper-pair size £ in uncon- 
ventional superconductors with a Russian-doll-like con- 
densate. For simplicity, assume that the Cooper-pair 
size is independent of temperature, or at lest, is a weak 
function of temperature. In real superconductors, if the 
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FIG. 8. Spatial variations of (a) the order parameter |^| 
and the averaged wavefunction \ip\ 2 , and (b) the energy gaps 
A p and A c inside and outside an isolated vortex in an infinite 
unconventional superconductor with a Russian-doll-like con- 
densate. R is the distance from the center of the vortex. £gl 
is the coherence length, and £ is the Cooper-pair size. For 
simplicity, all the functions are normalized at R — oo to 1. 
Since the pairing wavefunction ij) is not global, ij) in plot (a) 
is average. The letters A and B indicate a different behaviour 
of \4>\ 2 and Ap. 

pairing temperature T p is sufficiently high and there are 
structural transitions below T p , £ will vary with temper- 
ature. 

Let us consider the variations of the order parameter 
and the Cooper-pair wavefunction into the vortex cores 
appearing in the mixed state. The spatial variations of 
the order parameter | M/ 1 2 and the averaged wavefunction 
\vb\ 2 into a vortex core are sketched in Fig. 8(a). Inside 
a vortex core, the order parameter always goes to zero, 
while \ip\ 2 must not. Depending on the coupling strength 
of the Cooper pairs, there are two possible outcomes for 
the behaviour of \ip\ 2 into a vortex core, shown by the 
letters A and B in Fig. 8(a). The spatial variation of 
the averaged |i/>| 2 wavefunction, if such occurs, is of the 
order of £ in the clean limit. The coherence energy gap 
A c (cx 1^1) will be always zero inside a vortex core, as 
depicted in Fig. 8(b). In contrast, the pairing energy 
gap A p will not be zero into a vortex core, and even, 
may remain unchanged [the case A in Fig. 8(b)]. In 
the latter case, the size of the Cooper pairs cannot be 
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FIG. 9. Spatial variations of the total energy gap 
At = (Ap + Ac) 1 / 2 inside and outside an isolated vortex 
in an infinite unconventional superconductor with a Rus- 
sian-doll-like condensate (A c < A p ). The line R passes 
through the center of the vortex core. £gl is the coherence 
length, and £ is the size of Cooper pairs. The letters A and 
B correspond to the cases shown in Fig. 8(b). 

determined from the tunnelling data. In comparison with 
the case of two-band superconductors (Fig. 5), \i/j\ 2 and 
A p will never be zero inside a vortex: the presence of the 
long-range phase coherence in a superconductor with a 
Russian-doll-like condensate in the absence of the Cooper 
pairs has no sense. 

In superconductors in which the superconducting con- 
densate has a structure similar to that of a Russian doll, 

the total energy gap below T c is At = y^A 2 + A 2 [17]. 
Therefore, tunnelling measurements performed in vortex 
cores of such unconventional superconductors should the- 
oretically show the spatial variation of the total energy 
gap, as illustrated in Fig. 9. Since always A c < A p [17], 
the spatial variations of A t should be more pronounced 
within £ from the walls of a vortex (the cases B and C 
in Fig. 9). From the data, one can estimate not only the 
value of the coherence length but also the value of £. 

In real superconductors, however, the situation can 
be different. Let us consider an example: assume that 
in a superconductor with a Russian-doll-like condensate, 
magnetic fluctuations mediate the long-range phase co- 
herence. Tunnelling measurements probe the local den- 
sity of states of quasiparticle excitations present on the 
surface. It is known that, in magnetic materials, the 
spectrum of magnetic excitations on the surface differs 
from that in the bulk [4]. Depending on the quality of 
the surface, spin fluctuations into the topmost layer may 
not be able to mediate the long-range phase coherence. 
Therefore, the Cooper pairs on the surface will remain 
uncondensed. As a consequence, in tunnelling measure- 
ments performed in this superconductor, the pairing en- 
ergy gap Ap will be predominant in the spectra, and its 
spatial variations into a vortex core will be similar to 
those shown in Fig. 6(b). 
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B. Condensate with two types of Cooper pairs 

The second type of unconventional superconductors 
which we are going to discuss now is the superconduc- 
tors with a condensate consisting of two different types 
of Cooper pairs, which are independent of one another. 
These two types of Cooper pairs are formed due to dif- 
ferent bosonic excitations. The Cooper pairs of one type, 
for example, are formed due to phonons, while spin fluc- 
tuations, for instance, are responsible for the formation 
of Cooper pairs of the second type. The long-range phase 
coherence is mediated by the overlap of the wavefunctions 
of the Cooper pairs of one of these types, meaning that 
this type of the Cooper pairs appears at T c . Assume that 
the Cooper pairs of the other type are formed above the 
critical temperature, T p > T c . In these unconventional 
superconductors, the symmetries of the two types of the 
Cooper-pair wavefunctions can be different, as in the case 
of superconductors with a Russian-doll-like condensate. 

Let us consider two possible cases of the spatial dis- 
tribution of the Cooper pairs. If the two types of the 
Cooper pairs "penetrate" one another on a nanoscale, 
this situation is basically equivalent to the case of half- 
conventional superconductors, discussed above. There- 
fore, tunnelling spectroscopy data obtained along a line 
passing through the center of a vortex core in such super- 
conductors should be similar to those obtained in half- 
conventional superconductors. As a consequence, Fig- 
ures 4-6 are in principle applicable to this case, and the 
spatial variations of the two gaps are sketched in Fig. 
6. Depending on the quality of the surface, one can de- 
tect both gaps simultaneously or one at a time. From 
these data, one can in principle estimate the values of 
£ GL and/or £. 

If the two types of the Cooper pairs do not mix on a 
nanoscale, the superconducting phase will then have the 
"shape" of the Swiss cheese with the incoherent Cooper 
pairs of the other type inside the "cavities" (in layered su- 
perconductors, "pancakes"). On the surface, this struc- 
ture will look like the coexistence of two types of alter- 
nating patches. In the mixed state, the vortices will 
penetrate the superconductor in the weakest spots, i.e. 
in spots where the minimum free energy is needed to 
break the superconducting phase. In this case, the vor- 
tex lattice may be not regular. On the other hand, due to 
a strong interaction between vortices, the vortex lattice 
can become regular, and the "weakest superconducting 
phase" will follow the vortices. In any case, by perform- 
ing tunnelling measurements along a line passing through 
the center of a vortex core in such unconventional su- 
perconductors, one can get information only about one 
energy gap and, respectively, can estimate either £,gl or 
£2 (the size of the incoherent Cooper pairs). The spatial 
variations of the two energy gaps will be similar to those 
sketched in Fig. 6. 



V. VORTEX-CORE TUNNELLING DATA 
OBTAINED IN REAL SUPERCONDUCTORS 

In the literature, one can find only a few sets of tun- 
nelling spectroscopy data obtained in the mixed state 
of half-conventional and unconventional superconductors 
[5-14] . Here we discuss these data, expecting to estimate 
the values of the coherence length and the size of the 
Cooper pairs in respective superconductors. All super- 
conductors which will be considered in this section are 
obviously type-II superconductors. We start with half- 
conventional superconductors. All the data presented 
in this section are acquired in superconductor-insulator- 
normal metal junctions. 

A. Half-conventional superconductors 

The recently discovered magnesium diboride (MgB2) 
superconductor is a typical half-conventional supercon- 
ductor [4,17]. Ironically enough, it is probably the most 
studied half-conventional superconductor. 

In MgB2, superconductivity occurs in the boron layers. 
Band-structure calculations of MgB 2 show that there are 
at least two types of bands at the Fermi surface. The 
first one is a narrow band, built up of boron a orbitals, 
whilst the second one is a broader band with a smaller 
effective mass, built up mainly of it boron orbitals. 

The presence of two energy gaps in MgB2 is a well 
documented experimentally. The larger energy gap A CT 
occurs in the c-orbital band, the smaller gap A T in the 
7r-orbital band. The gap ratio 2A/ (fcsT c ) for A CT is about 
4.5. For Att, this ratio is around 1.7, so that, A^/A^ ~ 
2.7. Both energy gaps have an s-wave symmetry. The 
larger gap is highly anisotropic, while the smaller one 
is either isotropic or slightly anisotropic. The induced 
character of A T manifests itself in its temperature de- 
pendence: if A (j follows the temperature dependence de- 
rived in the framework of the BCS theory, the temper- 
ature dependence of A v lies below the BCS dependence 
at T -» T c . 

Tunnelling spectroscopy measurements have been per- 
formed in MgB2 in the mixed state along a line passing 
through the center of a vortex core [13]. Since the mea- 
surements have been conducted along the c-axis, they 
have been tuned outside of the vortex core to A„- which 
is the coherence gap in MgB2. As the authors admit, 
"the spectra in the center of a vortex are absolutely flat, 
..." [13]. This result is in good agreement with Fig. 6(a). 
The inner structure of a vortex would be complex if the 
measurements would be initially tuned to A CT . The coher- 
ence length £gl estimated from the size of vortex cores 
is about 500 A [13]. 

The behaviour of A^ in the presence of magnetic field 
is also documented in the literature. Figure 10 shows 
a set of tunnelling conductances recorded perpendicular 
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FIG. 10. Tunnelling conductances obtained in MgE>2 at 2 
K and different magnetic fields [14]. The spectra are recorded 
between the vortices. They are offset vertically for clarity. 
The critical temperature of the sample is 37.7 K. 

to the c-axis in applied magnetic fields up to 1 T. The 
conductances are obtained between vortices, not inside 
of a vortex. In Fig. 10, the conductance at zero field has 
a double-gap structure reflecting the presence of both 
A,,- and A^. As the field is increased, the inner peaks 
associated with A T are rapidly suppressed, and the zero- 
bias conductance increases. The outer peaks reflecting 
A a are affected by the field much less than those from 
A^. In the absence of the field, A CT ~ 7.2 meV, and at H 
= 1 T A CT ~ 6.6 meV. If these spectra were taken along 
a line passing through the center of a vortex, they would 
correspond to the case B in Fig. 6(b). 

The second half-conventional superconductor which 
was also intensively studied in the mixed state by tun- 
nelling spectroscopy is 2H-NbSe2 [5-7]. 2H-NbSe2 is 
a "clean" layered superconductor with the presence of 
charge-density waves (thus, it is low-dimensional). Be- 
low T c , 2H-NbSe2 has two energy gaps, and the ratio 
A^/As is about 3 [23,24]. As a consequence, £gl/£ ~ 3 
too. The conductances obtained in the center of a vortex 
exhibit a a peak at zero bias which will be discussed be- 
low. The coherence length ^gl estimated from the size 
of vortex cores is about 175 A. 

B. Unconventional superconductors 

In this subsection, we shall discuss tunnelling spec- 
troscopy data obtained in three unconventional supercon- 
ductors. We start with the data measured in LUNI2B2C 
[9]. 

The nickel borocarbide class of superconductors has 
the general formula i?Ni2B2C, where R is a rare earth 
which is either magnetic (Tm, Er, Ho, or Dy) or nonmag- 
netic (Lu and Y). The Ni borocarbides have a layered- 
tetragonal structure alternating RC sheets and Ni 2 B 2 
layers. Transition temperatures in these quaternary in- 




Bias (mV) 

FIG. 11. Tunnelling conductances obtained in LuNi2B2C 
at 4.2 K and 0.375 T [9]. The spectra are recorded at 470 
A (curve A), 130 A (B) and 20 A (C) from the center of a 
vortex. The curve D is measured in the center of the vortex. 
The critical temperature of the sample is 15.8 K. 

termetallic compounds can be as high as 17 K. In the 
case when R = Pr, Nd, Sm, Gd or Tb in i?Ni 2 B 2 C, the 
Ni borocarbides are not superconducting at low tempera- 
tures but antiferromagnetic. In the Ni borocarbides with 
a magnetic rare earth, superconductivity coexists at low 
temperatures with a long-range antiferromagnetic order. 

Many different types of measurements carried out in 
the Ni borocarbides show that the gap ratio 2A/(&bT c ) 
is between 3.2 and 5.3 [17]. What concerns the shape of 
the energy gap, there is complete disagreement in the lit- 
erature. In photocmission and microwave measurements, 
the energy gap in some Ni borocarbides was found to 
be an s-wave but highly anisotropic [17]. On the other 
hand, in specific-heat, thermal-conductivity and Raman- 
scattering measurements carried out in the Ni borocar- 
bides with R = Y and Lu, the energy gap was found to 
be a highly anisotropic gap, most likely with nodes. Fur- 
thermore, in other thermal-conductivity measurements, 
the gap appears to have point nodes along the [100] and 
[010] directions, thus along the a and b axes. Recent tun- 
nelling measurements performed in the antiferromagnetic 
TmNi2B2C unambiguously show that this Ni borocarbide 
is a fully gapped s-wave superconductor with a gap be- 
ing slightly anisotropic. To reconcile all these data, one 
should assume that different measurements probe differ- 
ent energy gaps, either A p or A c . 

Figure 11 depicts tunnelling conductances obtained in 
the mixed state in the LuNi2B2C borocarbide at various 
distances from the center of a vortex. All the conduc- 
tance in Fig. 11, including that recorded in the center of 
a vortex, have a gap structure. However, the humps in 
the conductances obtained inside the vortex are not sit- 
uated symmetrically relatively zero. The bias positions 
of the left-hand humps in the spectra are VJ(470 A) ~ 
3.4 mV, V;(130 A) ~ 3.1 mV, V t (20 A) ~ 2.5 mV, and 
Vi(0 A) ~ 1.8 mV. The bias positions of the right-hand 
humps in the spectra are V r (470 A) ~ 3.4 mV, K(130 A) 
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FIG. 12. The phase diagram of hole-doped cuprates 
[25,26,4,17]. A c is the phase-coherence energy gap, and A p is 
the pairing energy gap. 

~ 2.9 mV, U r (20 A) ~ 3.0 mV, and V r (0 A) ~ 3.1 mV. 
The bias positions of the humps are determined with an 
error of about ±0.1 mV. Unfortunately, from these data, 
it is impossible to determine what type of unconventional 
superconductors LuNi 2 B 2 C belongs to (see the previous 
section). The evolution of the hump position on the left- 
hand side of the conductances is in agreement with the 
case A in Fig. 9. At the same time, the evolution of the 
hump position on the right-hand side of the conductances 
is in agreement with the case B in Fig. 6(b). Indeed, the 
gap ratio 2A(470 A)/fcsT c ~ 4.9 indicates that the spec- 
tra in Fig. 11 reflect either a A p (A^) or a A t gap (see 
above). The radius of the vortex cores in LuNi 2 B 2 C is 
about 200 ± 10 A. Then, £ Gi « 200 A. If we rely on 
the bias position of the right-hand humps, then one can 
estimate the size of the Cooper pairs: £ ~ 200 - 130 = 70 
A. If this is the case, then A^/As oc £gl/(, ~ 3. Unfor- 
tunately, the set of the data presented in Fig. 11 is too 
limited for more definitive evaluation of the size of the 
Cooper pairs. 

Let us now discuss tunnelling data obtained in the 
mixed state of two cuprates: YBa2Cu307_«5 (YBCO) and 
Bi 2 Sr 2 CaCu208 (Bi2212). The superconducting phase in 
cuprates appears when they are slightly doped. The crit- 
ical temperature of cuprates can be tuned by varying the 
doping level. Hole-doped cuprates have two energy gaps 
related to the superconducting phase, A p and A c , shown 
in Fig. 12. The electron-doped cuprates have a similar 
phase diagram [27,4]. 

Figure 13 shows two averaged tunnelling conductances 
recorded in the mixed state of Bi2212. The conductance 
obtained in the center of a vortex has two features: an 
apparent enlargement of the gap magnitude and the pres- 
ence of a subgap. The subgap corresponds to the gap 
inside the cortex cores, the magnitude of which is pro- 
portional to the magnitude of the gap recorded outside of 
the vortices with the coefficient of about 0.3 [12]. Thus, 
the subgap in the lower conductance in Fig. 13 is either 




FIG. 13. Averaged spectra recorded in the center of a vor- 
tex and outside the vortex in overdoped Bi2212 (T c — 87.4 K) 
at 4.2 K and 6 T [12]. The arrows indicate a subgap structure. 
The upper curve is offset vertically for clarity. 

A p [the case B in Fig. 6(b)] or A t [the case B in Fig. 9] 
with a reduced magnitude. From the data presented in 
[12], it is very difficult to distinguish between these two 
possible cases: seemingly, the value of the subgap remain 
constant near to the center of the vortex (see Fig. 2(a) 
in [12]). So, it is more likely that the subgap in the lower 
conductance in Fig. 13 corresponds to A p . The radius 
of the vortex core in slightly overdoped Bi2212 is about 
22 ± 3 A [11]. Because of the small size of the vortices 
in Bi2212 and a weak manifestation of the subgap in the 
spectra [12], it is practically impossible to estimate the 
size of the Cooper pairs in Bi2212 from the data. The 
only way to estimate £ is from the ratio £/£gl ~ A c /A p . 
Then, by using Fig. 12 one can obtain that in slightly 
overdoped Bi2212, £ - 22 x 20/32 A ~ 14 A. 

The increase of the magnitude of the gap in Fig. 13 is 
apparent. The quasiparticle peaks in tunnelling conduc- 
tances recorded below T c in cuprates appear on top of 
a contribution from a normal-state pseudogap [28,29,4]. 
The disappearance of the peaks discloses the humps from 
the pseudogap. Since the magnitude of the pseudogap is 
larger than the magnitudes of A p and A c , this creates an 
effect of an increase of the magnitude of the gap inside 
the vortex cores. This can be demonstrated by using the 
data from [10]: in Fig. 14, one can see that the conduc- 
tance recorded in the center of a vortex is very similar to 
the conductance obtained above T c . The gap structure in 
both these curves is caused by a normal-state pseudogap 
which is always present "beneath" the superconducting 
gap(s). The same effect of an apparent increase of the 
magnitude of the gap inside the vortex cores also occurs 
in YBCO, as illustrated in Fig. 15. 

The subgap structures in conductances recorded in the 
center of a vortex are more pronounced in YBCO than 
those in Bi2212, as shown in Fig. 15. The problem, 
however, is that this subgap is caused not only by the 
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FIG. 14. Tunnelling conductances obtained in underdoped 
Bi2212 with T c = 83 K [10]. The lower conductance is taken 
in zero magnetic field at 98.6 K. The other curves are mea- 
sured at 4.2 K and 6 T. The upper curve is recorded between 
vortices, and the middle curve is acquired in the center of a 
vortex. For clarity the 4.2 K curves are offset vertically. 
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FIG. 15. Tunnelling conductances obtained in an YBCO 
single crystal with T c — 91 K at 4.2 K [8]. The upper curve 
was measured in zero magnetic field. The lower curve is the 
averaged spectrum of conductances obtained at the center of 
a vortex core along a 5 nm path and in H — 6 T applied 
along the c axis. The lower curve was magnified by a factor 
of 1.875. The upper curve is offset vertically for clarity. 



superconducting gap with a reduced magnitude, as that 
in Bi2212, but also by the gap on chains. In the crys- 
tal structure, the unit cell of YBCO has four CuO chains 
which are parallel to the b crystal axis. YBCO is the only 
superconducting cuprate having one-dimensional CuO 
chains. In YBCO, the CuO chains become supercon- 
ducting due to the proximity effect [17]. The value of 
the superconducting energy gap on the chains in opti- 
mally doped YBCO is about 6 meV [4]. In Fig. 15, one 
can see that even the conductance obtained without an 
applied magnetic field has a subgap structure caused by 
the superconducting energy gap on the chains. In the 
vortex cores, the bulk superconducting gap with a re- 
duced magnitude appears in tunnelling spectra "on top" 
of the superconducting gap on CuO chains. Because of 
this, it is impossible to estimate of the Cooper-pair size 
from the data presented in [8]: the weak subgap in con- 
ductances recorded outside a vortex smoothly transforms 
to a more-pronounced subgap in conductances acquired 
in the vortex core, having a very similar magnitude (see 
Fig. 3(b) in [8]). In YBCO, the size of the vortex is at 
least a factor of 2 larger than that in Bi2212 [10], thus, it 
is slightly larger than 40 A. However, the vortex size in 
the bulk of YBCO, determined by /iSR at 6 T, is about 
20 A [3]. It is also worth to mention that, in Fig. 15, the 
magnitude of the energy gap in the conductance recorded 
without an applied magnetic field corresponds to A c in 
Fig. 12, and not to A p , as those in Figs. 13 and 14 in 
Bi2212. 



C. Zero-bias conductance peak 

Tunnelling conductances obtained in the center of a 
vortex in 2H-NbSe2 exhibit a peak at zero bias [5-7] 
shown in Fig. 16. It is not the purpose of this paper 
to discuss the origin of this peak. However, it is worth to 
spent some time on this issue. Generally speaking, the 
zero-bias conductance peak is a manifestation of some 
bound states inside the vortex cores. The question is 
what kind of bound states are they? 

As shown in [4] , solitonic states present in a system will 
manifest themselves in tunnelling conductances through 
the appearance of a peak at zero bias. In superconductor- 
insulator-normal metal junctions, the shape of this zero- 
bias conductance peak will be proportional to a sech 2 
function as follows 

^l=Axsech 2 (V/V ), (5) 

where V is the applied bias, and A and Vb are the con- 
stants. The corresponding I(V) characteristic is 

I{V) = AV x tanh(VyFo). (6) 

It is interesting that the zero-bias conductance peak 
recorded in the center of vortex cores in 2H-NbSe2 also 
has the shape of a sech 2 function, as depicted in Fig. 
16. Superconductivity in 2H-NbSe2 coexist with charge- 
density waves which are quasi-one dimensional. Solitons 
require the presence of quasi-one-dimensionality in the 
system. So, it is possible that the zero-bias peak in con- 
ductances recorded in the center of vortex cores in 2H- 
NbSe2 is a manifestation of solitonic states in 2H-NbSe 2 
(which may be paired outside the vortex cores). 
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FIG. 16. Tunnelling conductance recorded in the center of 
a vortex in 2H-NbSe 2 and the v4xsech 2 (V/Vb) fit. The data 
are taken from Fig. 2 in [6]. In the fit, A = 2.5 and V = 1/2.8 
mV. 



A peak at zero bias also appear in conductances 
recorded in non-superconducting and superconducting 
compounds without an applied magnetic field. In the 
charge-density- wave conductor NbSe3, a zero-bias con- 
ductance peak has the shape of a sech 2 function too 
[30]. In superconducting hole-doped [4] and electron- 
doped [27,4] cuprates, a peak at zero bias occurring in 
conductances recorded usually at 45° relative to the in- 
plane crystal axes has also the shape of a sech 2 function. 



VI. DISCUSSION 

It is necessary to note that the determination of the 
values of the coherence length and the size of the Cooper 
pairs in unconventional superconductors by using tun- 
nelling spectroscopy is not a straight forward procedure. 
In unconventional superconductors, the values of the co- 
herence length and the size of the Cooper pairs are usu- 
ally small. As a consequence, the values of £gl and, in 
particular, £ estimated from the experiment will have a 
large error. In addition to this, there are at least three 
other aspects affecting the results. First, depending on 
the crystal structure, the shape of vortex cores can be not 
round: in 2H-NbSe2, the vortex cores have the shape of a 
six- leg star [6]. Second, the size of vortex cores varies: it 
depends on the distance between the vortices, i.e. finally, 
on H. Third, the vortices slightly expand near the sur- 
face. The size of a vortex can increase by up to 30% at 
the sample surface in comparison with that in the bulk 
[3] . A discussion concerning the structure and the size of 
vortex cores can be found elsewhere [31,3]. 

In the framework of the Ginzburg-Landau theory, the 
value of the coherence length in type-II superconductors 
can be determined from iJ c2 as follows 



$ 



2ttH, 



c2 



(7) 



where 



h 

~2e 



2.0679 x 10~ 15 Tm 2 (or Weber) 



(8) 



is the magnetic flux quantum. In layered unconventional 
superconductors, the critical magnetic field H C 2, as well 
as £gl, is different in different directions — parallel and 
perpendicular to the layers. For example, the upper 
critical field applied perpendicular to the layers, H c2y ±, 
is determined by vortices whose screening currents flow 
parallel to the planes. For layered superconductors, the 
Ginzburg-Landau expression takes the following form 



H C 2.± — 



4> n 



2*e 



(9) 



GL.ab 



where the letters "ab" indicate that the direction of the 
screening currents is in the ab-plane. The value of the 
coherence length perpendicular to the planes, £gl,c can 
be determined from the ratio 



H C 2,\\ _ S,GL,ab 
H C 2,± ^GL.c 



(10) 



where H c2 ^\ is the upper critical field applied parallel to 
the layers. 

In reality, however, the values of £,GL,ab and £gl,c> 
determined in unconventional superconductors through 
H c 2 t \\ and if C 2,J- are not the in-plane and out-of-plane 
coherence lengths — they correspond to the size of the 
Cooper pairs in two directions, i.e. to £ a f, and £ c [4]. 
As an example, let us consider the doping dependence 
of the in-plane coherence length and Cooper-pair size in 
La2- 2; Sr 2 ;Cu04 (LSCO). The lower curve in Fig. 17 is 
the in-plane characteristic length in LSCO determined 
from the -ff C 2,_i_(p) dependence. It is easy to check that 
the dependence l/£ a b(f>) from Fig. 17 does not mimic 
the A c (p) dependence in Fig. 12. Instead, it is similar 
to the A p (p) dependence in Fig. 12. In fact, the up- 
per curve in Fig. 17 represents £gl(p) in LSCO which 
will be discussed in a moment. Let us first consider 
other examples: in MgB2, the value of the coherence 
length determined from the size of vortex cores is about 
500 A (see above), while determined from ff c2 is about 
100 A [13]. In 2H-NbSe 2 , £gl ~ 175 A, while £ ~ 
77 A, respectively [6]. In LuNi 2 B 2 C, £ GL ~ 200 A, 
while £ ~ 75 A, respectively [9]. In slightly overdoped 
Bi2212, £ GL ~ 22 A (see above), while £ ~ 14 A de- 
termined from H C 2 — $o/(27r£ 2 ). All these values are 
low-temperature values. So, experimentally, the value of 
the coherence length estimated from the vortex-core size 
is always larger than that obtained from the Ginzburg- 
Landau expression H c2 = $0/ '(-^Cgl)- ^ n unconven- 
tional superconductors, the H C 2 method provides the val- 
ues of the Cooper-pair size. 

In unconventional superconductors, it is possible to es- 
timate the values of £ and £gl independently. For ex- 
ample, the Cooper-pair size at low temperature can be 
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FIG. 17. Sketch of doping dependences of low-temperature 
in-plane £gl(0) and £(0) in LSCO. The Cooper-pair size is 
obtained from H c2 = $ /(2tt^ 2 ) [32]. 

obtained from Eq. (1). To estimate the size of vortex 
cores is more complicated — such an approach is described 
elsewhere [32]. Figure 17 depicts the doping dependence 
of ^GL,ab in LSCO, derived by using the doping depen- 
dences of other superconducting characteristics, such as 
the critical current density, the collective pinning energy, 
the supernuid density and the condensation energy. In- 
deed, the l/£,GL,ab(p) dependence from Fig. 17 is similar 
to the doping dependence A c (p) in Fig. 12. At p ~ 1/8, 
£,GL,ab(p) has a kink related to the so-called g anomaly 
inherent exclusively to LSCO. This fact indicates that the 
disappearance of superconductivity in LSCO at p ~ 1/8 
is due to the absence of the long-phase coherence, not 
the Cooper pairs. 

VII. CONCLUSIONS 

The main purpose of the paper was to discuss a possi- 
bility of the determination of the values of the coherence 
length and the Cooper-pair size in unconventional super- 
conductors by using tunnelling spectroscopy. In uncon- 
ventional superconductors, the inner structure of a vortex 
core has a complex structure determined by the order pa- 
rameter of the superconducting state and by the pairing 
wavefunction of the Cooper pairs. In clean superconduc- 
tors, the spatial variations of the order parameter and 
the pairing wavefunction occur over the distances of the 
order of the coherence length and the Cooper-pair size, 
respectively. Therefore, by performing tunnelling spec- 
troscopy along a line passing through a vortex core, one 
is able, in principle, to estimate the values of the coherent 
length and the Cooper-pair size. In the paper, common 
patterns of the structure of the vortex cores have been ob- 
tained for two-band (half-conventional) superconductors 
and unconventional superconductors. At the moment, a 
detail comparison tunnelling spectroscopy data recorded 
in half-conventional and unconventional superconductors 



and the derived patterns is not possible because the set of 
tunnelling spectroscopy data available in the literature is 
limited. In the future, both the spatial and energy reso- 
lutions of tunnelling spectroscopy measurements need to 
be improved since the size of the vortex cores in uncon- 
ventional superconductors is very small. 
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